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In this paper, we present an approximate expression for determining the effective permittivity 
describing the coherent propagation of an electromagnetic wave in random media. Under the Qua- 
sicrystalline Coherent Potential Approximation (QC-CPA), it is known that multiple scattering 
theory provided an expression for this effective permittivity. The numerical evaluation of this one 
is, however, a challenging problem. To find a tractable expression, we add some new approximations 
to the (QC-CPA) approach. As a result, we obtained an expression for the effective permittivity 
which contained at the same time the Maxwell-Garnett formula in the low frequency limit, and the 
Keller formula, which has been recently proved to be in good agreement for particles exceeding the 
wavelength. 

PACS numbers: 42.25.Bs, 41.20.Jb, 78.20.-e 



I. INTRODUCTION 

The description of electromagnetic waves propagation 
in random media in term of the properties of the con- 
stituents has been studied extensively in t he p ast decades 
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llllll|2fl|2l|2l|2l|2||- In most of works, the ba- 
sic idea is to calculate several statistical moments of the 
electromagnetic field to understand how the wave inter- 
act with the the random medium |^ ITTl IT3 . IT3 . H^ . 
In this paper, we are concerned by the first moment 
which is the average electric field. Under some assump- 
tion, it can be shown that the average electric field 
propagates as if the medium where homogeneous but 
with a renormalized permittivity, termed effective per- 
mittivity. The calculation of this parameter as a long 
history which dates back from the work of Clausius- 
Mossotti and Maxwell Garnett f2^. Since then, most 
of the study are concerned with the quasi-static limit 
where retardation effect are iieglected |2ll 12^. l2^ l2^ 
lilllllillllllillilsllllllllll. In order to take 
into account scattering effects, quantum multiple scat- 
tering theory has been transposed in the electromagnetic 
case 0, 0, 0, 0, 01 , but as a rigorous analytical an- 
swer is unreachable^several approximation schemes have 




is the Quasicrystalline Coherent Potential Approxima- 
tion (QC-CPA) which takes into account the correlation 
between the particles |il3, J6, ,3^ J£, ^0] . Un- 

fortunately, except at low frequency, the answer is still 
too involved to permit the calculation of the effective 
permittivity. 

The aim of this paper is to add some new approx- 



imations to the (QC-CPA) approach which furnish a 
tractable equation for the effective permittivity. The ex- 
pression obtained contains the low frequency limit of the 
(QC-CPA) approach. At this limit, the (QC-CPA) equa- 
tions can be written as a generalized Maxwell Garnett 
formula and are proven to be in g ood agreement with the 
experimental results 0, 0, 113, ■ Furthermore, 
the formula obtained contains also the approximate for- 
mula due to Keller, which has been derived in using scalar 
theory, but seems to be in accord with the experimental 
data iFor particles larger than a wavelength j52. IssL Is^ . 

The paper is organized as follows. In Sectional we 
introduce the multiple scattering formalism and we show 
under what hypothesis the effective medium theory is 
valid. In section lTTll we recall the different steps in order 
to obtain the system of equation verified by the effec- 
tive permittivity under the (QC-CPA) approach. Then, 
we introduce, in section Hvl some new approximations in 
order to obtain a tractable formula for the effective per- 
mittivity. In the two following section IVII we derive 
respectively the low frequency and high frequency limit 
of our new approach. 



II. DYSON EQUATION AND EFFECTIVE 
PERMITTIVITY 



In the following, we consider harmonic waves with 
e~'"* pulsation. We consider an ensemble of 3> 1 iden- 
tical spheres of radius with dielectric function £^(0;) 
within a infinite medium with dielectric function (.i{uj). 
The field produced at r by a discrete source located at Tq 
is given by the dyadic Green function G(r, Tq, lS), which 
verifies the following propagation equation: 

V X V X G(r,ro,w) - ev{r,uj) Kl^^G{r,n^,uj) 

= 6(r - rn)I (1) 
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where Kyac = t^/c with c the speed of hght in vaccum 
and 

N 

ey(r,cj) = ei(w) +^[es(w) - ei(w)] O^Cr - Tj) , 
i=i 

where ri, . . . ,r]\[ are the center of the particles and Qg 
describes the spherical particle shape : 



e.(r) 



1 if ||r|| < rs 
if llrll > r. 



(2) 



The solution of equation |^ is uniquely defined if we 
impose the radiation condition at infinity. 

The multiple scattering process by the particles is 
mathematically decomposed in introducing the Green 
function Gi , describing the propagation within an ho- 
mogenous medium with permittivity £1(0;), which verifies 
the following equation: 

V X V X G^{r,ro,uj) - ei{uj) K^^^G'^ {r,ro,uj) 

= 5{r -ro)7, (3) 

with the appropriate boundary conditions. In an infinite 
random medium, we have |l6l Issl Is^ : 



Gi ir,ro,Lu) = 



- VV 



iKi 11? — rn\\ 

47r||r|| 



(4) 



where = £1(^)^-2^^. 

In using this Green function, we decompose the Green 
function G{r,ro,ui) under the following form |^ ITlI 

M 

G = gT + gT V G, (5) 
where the following operator notation is used: 

[A-B]{r,ro)^ / d^ n A(r, n) • B(ri, ro) . (6) 



The potential V, which describes the interaction between 
the wave and the particles, is given by: 



N 



(7) 



(r, ro, = (27r)2 S{r-ro) f r. {r, uj) , (8) 
Vr, {r, ij) = [i^f - Kf] Bdir - n)I . (9) 

with Kg — es(uj) K'i„^. It is useful to introduce the T 
matrix defined by jSllllIllIl: 



g = gT + gT-t.gI 



(10) 



In iterating equation ^ and comparing it with the defini- 
tion 11U|) , we show that the T matrix verifies the following 
equation: 



T + V G^ T . 



(11) 



If we introduce the T matrix of each scatterer by: 



(12) 
(13) 



we can decompose the T matrix for the whole system, in a 
series of multiple scattering processes by the particles 

millLie]: 



N N N 

T = ^ ^ tj.^ -\- ^ ^ ^ ^ trj • G^ • tj' 



(14) 



This T matrix is useful to calculate the average field < 
G > since we have: 



< G >= Gi +Gi • < T > -Gi 



(15) 



The equivalent of the potential operator V for the aver- 
age Green function < G > is the mass operator S defined 
by: 



< G >= Gi + Gi ■ S- < G > 



(16) 



Similarly to equation Hll|) . we have the following rela- 
tionship between the average T matrix and the mass op- 
erator: 



<r>=S + S-Gi ■ <T> 
or equivalently, 



S =< T > • 



/ + Gi < T > 



(17) 



(18) 



The mass operator correspond to all irreducible diagrams 
in the Feynman representation |3, ^| . The equa- 
tion (|16|) written in differential form is: 

V X Vx <G{r,ro,uj) > ~ei{uj) KI^^ <G{r,ro,uj) > 

- J d^ri S(r,ri,u;)- < G{ri,rQ,uj) >= 6{r - vq)! . 

(19) 

For a statistical homogeneous medium we have: 

S(r,ri,cj) = S(r-ri,w), (20) 
< G(r, ro, Lo) > G(r - ro, lo) > . (21) 

Thus, we can use a Fourier transform: 

S(fc,w) = J d^r exp(-ifc • r) S(r,w) , 
G{k,uj) = J d^r exp(-i fc-r)G(r,w), 

and equation (|19l) becomes: 

k\\^ (I - kk) ~ ei{Lu) K^^J -^{k,uj)] ■ <G{k,uj) >=J. 

(24) 



(22) 
(23) 
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For a statistical isotropic medium, we have: 

S(fc,tj) = Si(||fc||,cj)(7-fcfc) + S||(||fc||,cj)fcfc. (25) 
with k = fc/||fc|| and then: 



Furthermore, we see that the contribution due to pole 
i^T = in the second term of equation (|32|l is null; in 
fact, the dyadic VV operate on a constant since we have 
giif ||r|| _ 2 fQj. ^]-^jg pole. Hence, we obtain the following 
expression for the Green function: 



< G(fc, Lo)>= \\k\mi- kk) - ei (lu) Klj 



I kk 



S(fe,w) 

(26) 



fcfc 



S±(||fc||,^)) 



£1 



(^)if2^, + I]||(||fe||,^) 



(27) 



In the following, we introduce two effective permittiv- 
ity function and el defined by: 



tl{\\k\Uuj)Kl^^ = ei{uj)Kl^^- 
and H27|l is written: 
<G{k,u) > 

kk 



\k\luj) 



(28) 
(29) 



kk 



\k\\^-eHmiio)K' 
kk 



2 

vac 



4i\M,uj)K^,, el{\\k\\,to)K^,, ' 
The Green function in the space domain is: 



(30) 



<G{r,Lj) > 
d^k 
(2^ 

d^fc 

(27r)3 



VV 



kr 



i\\k\\,uj)K^^ 



kk 



4mi^)KLc el(\\k\luj)KL 



|fc||2 • 
(31) 



After an integration on the solid angle in equation l|31() 
given by the expression ljAl|l in the appendix, we obtain: 



< G{r,Lu) > 



1 



i\\r\ 
1 

illrl 



+ 00 



dK 
(2^ 



VV 



< G{r,u) >-- 



E 

1=1 



- VV 



.i^e. \\r\\ 

47r||r|| 



(33) 



where Kei are the roots of K^^ = e-^{Kei,u})K^ac such 
as lm{Kei) > to insure that the radiation condition 
at infinity is verified. Sheng has called the roots Kei 
the quasi- modes of the random media 0, 0. If we 
only consider the root Kg = Kej which has the smallest 
imaginary part {lu\{K(.j) = mini [Ini(^ei)]) a-nd then 
the smallest exponential factor in equation 133|) . we de- 
fine the effective permittivity by ee(w) = e^(ife,w). The 
average Green function is then equal to the Green func- 
tion for an infinite homogenous medium with permittiv- 
ity ee{uj) : 



< G{r,u) >=G^(r,w), 



where 



VV 



e^K, \\r\\ 

47r||r|| 



(34) 



(35) 



Thus, the effective medium approach is valid if we ne- 
glect the longitudinal excitation in the medium and if 
the propagative mode with the smallest imaginary part 
is the primary contribution in the developpement H33|) . 



III. THE COHERENT-POTENTIAL AND 
QUASI-CRYSTALLINE APPROXIMATIONS 

Previously, we have shown how the mass operator is re- 
lated to the effective permittivity. To calculate the mass 
operator, we can use equations ((njl and l(TS|) . However, 
we can improve this system of equations in rewriting the 
Green function development |SJl in replacing the Green 
function G^ by Gg : 



G = Gg +Gg Ve G, 
where we have to introduce a new potential V e'- 



-e^{K,uj)K, 



2 

vac 



N 



rVV 



dK 



iK\\r\\ 



\K\\r\ 



ei{K,u)K^,, el{K,Lo)K^^^ 

(32) 



1 

K 



i=l 



where we have supposed that e^(||fe||,w) = e^(-||fc||, cj) 
and e|!(||fc||,ti;) — ee(-||fc||,u;) . In using the residue the- 
orem, we easily evaluate these integrals. However, we 
neglect the longitudinal excitation, which are solutions 
of e-^{K,uj) = and el{K,uj) = 0, since we are only 
interested by the propagation of the transversal field. 



Vc,r, {r, To, w) = (27r)^ 5{r - Tq) ■Ue,r, (r, uj) , 
Ve^,^ (r, uj) - [Kl - Kl] Qs{r - r,)I 



Kl]I. 



(36) 

(37) 
(38) 

(39) 



Similarly to the previous section, we introduce a T matrix 
such that: 



G = G. 



G„ ■ T c ■ G„ 



(40) 
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which admits the following decomposition: 

N N N 



g: 



(41) 



where we have defined a renormalized T matrix for the 
particles: 



g: 



(42) 



In supposing that the effective medium approach is cor- 
rect, we impose the following condition on the average 
field: 



< G{r,uj) >== Ge(r,tj) , 



or equivalently, 



(43) 



(44) 



due to equation H4U|) . The condition 1438 is the Coherent- 
Potential Approximation (CPA) [H^lllllllli^. The 
expression (Hit and lUT)) form a closed system of equa- 
tions on the unknown permittivity ee(w). To the first or- 
der in density of particles, this system of equations gives 
equation: 



N 
i=l 



>= 0. 



(45) 



In Fourier-space, the T matrix for one scatterer verifies 
the property: 



te,r,(fc|fco) =e-'('=-'=«)-'-' *e,o(fc|fco), 



(46) 



where te,o is the T matrix for a particle located at the ori- 
gin of coordinate. The average of the exponential term, 
introduced in equation (|45|l by the properties H46|) . gives 
for a statistical homogeneous medium: 



N 

1=1 



N 



d3T._e-i(fe-fco) r ^ (47) 



n (27r)2 (5(fc - fco) 



(48) 



where we have introduced the density of scatterers n = 
N/V with V the volume of the random medium. The 
condition (|45() becomes: 



*e,o(fco|feo) = 0. 



(49) 



This CPA condition has been used in several works 0, 
0,|53|. It is worth mentioning that operator te^riik\ko) 
is not the T matrix describing the scattering by a particle 
of permittivity es{uj) surrounded by a medium of permit- 
tivity ee{uj). To do so, the operator (|39|l should have the 
following form: 



(50) 



However, we see that the operator Ve.ri {r, w) is quiet 
different from the equation H39(l . and in particular we 



have Ve-Ti (^j '-^ 



[Kf — Kl]I for I |r — 1 1 > contrary 
to the definition H5UI) . where Ve;^. (r, (jj) — when r is 
outside the particle. Thus, the operator t^^n is non-local 
and cannot be obtained from the classical Mie theory 0, 
l^iH^- To overcome this difficulty, the operator t^^n is 
replaced by the scattering operator of a " structural unit" 
in the works [Til I21L ISTj . Nevertheless, this approach 
doesn't seem to have any theoretical justification. 

Hence, we prefer to use the more rigorous approach 
introduce in the scattering theory by disorder liquid 
metal [s^ [s^ and adapted in the electromagnetic case 
by Tsang et al. [Hill. In this approach, the non-local 
term [Kf — K^]! is correctly taking into account by av- 
eraging equations (|41|) where the correct potential v^^n , 
defined by H39() . is used. A system of hierarchic equa- 
tions is obtained where correlation functions between two 
or more particles are successively introduced. The chain 
of equations is closed in using the Quasi- Crystalline Ap- 
proximation (QCA), which neglect the fluctuation of the 
effective field, acting on a particle located at rj, due to a 
deviation of a particle located at r; from its average po- 
sition [36|. This approximation describes the correlation 
between the particles, only with a two-point correlation 



function g{ri, rj ) 



9{ 



Under the QC-CPA 



scheme, we obtain the following expression for the mass 
operator [Illillll 38, 39]: 



S(fco,t^) = nCe,o{ko\ko) , 
Ce,o{k\ko) = te^o{k\ko) 
d^fci 



(51) 



(27r)3 



hik - fei) t,,o{k\ki) ■ G, (fci) • Ce,o(fcl|feo) 

(52) 



where 



• G„ ■ tf 



VeA'^^'^o) = {2TT)S{r ~ ro)tTe,o(r) , 



(53) 
(54) 
(55) 



and 



h(r) = g{r) - 1 , (56) 
h{k -ki)= j d^r exp(-i (fc - fci) ■ r) h{r) , (57) 



G7{k) 



d"^r exp(— ifc • r) G^ (r) . 



(58) 



If we rewrite the potential H55|l under the following form: 



(59) 



where we have defined a new wave number = — 
Kf + K'^, we see that the operator t^ o is the T matrix for 
a scatterer of permittivity = — ei + in a medium 
of permittivity eg- 
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As it is described in the previous section, the effec- 
tive propagation K,. constant is the root, which has the 
smallest imaginary part, of the equation: 

K^^K^ + ^^{K,,io), (60) 

where the mass operator is decomposed under the form 
(|25|l . Once the effective wave number obtained, the 
effective permittivity is given by: 

eeH^i^eV^ic- (61) 

IV. SOME FURTHER APPROXIMATIONS 

As it can be guessed, solving numerically the pre- 
vious system of equations (|51l61f) is full of complex- 
ities. However, the low frequency limit of this sys- 
tem of equation has been obtained analytically and has 
shown to be in good agreement with the experimental 
results 13, 16]. We have also to mention that the numer- 
ical solution of the quasicrystalline approximation (but 
without the coherent potential approximation) has been 
developed [HiOI. 

To reduce the numerical difficulties in the system of 
equations (|51I61|) . we add two new approximations to 
the QC-CPA scheme: 

• A far-field approximation: For an incident plane 
wave: 

E\r) = E'{ko)e' , (62) 

transverse to the propagation direction fco^ 

E\ko)-ko = 0, (63) 

where ko = Keko and ko ■ ko = 1, the scattered 
far-field, by a particle within a medium of permit- 
tivity ee((x!), is described by an operator /(fc|fco) 
such that: 

E%r) = ^^^l^7{k\ko) ■ E\ko) . (64) 

which verifies transversality conditions: 

7(fe|fco)-fco = 0, (65) 
fc-7(fc|fco)=0. (66) 

Moreover, the scattered field in the general case is 
expressed with the operator te^o by: 

£;^(r)= J d^nd^r^G^ {r,r,) ■t,,,o{ri\r2) ■ E\r2) . 

(67) 

In using the phase perturbation method in equation 
l)67|l . the scattered far-field is obtained in function 
of the operator te,o{k\kQ), and in comparing the 



result with equation H64|l . we obtain the following 
relationship: : 

47r7(fc|feo) ={I-kk) 

■te.o{Kek\Keko)-(l-koko), (68) 

Our far-field approximation consist in neglecting 
the longitudinal component and the off-shell con- 
tribution in the operator tg.o, and we write: 

t(i^efc|ifefco) -4^7(fe|fco), (69) 

= 4^ (I - kk) -Jiklko) ■ (I - koko) ■ 

(70) 

where the last equality comes from the properties 

• A forward scattering approximation: For scatter- 
ers large compared to a wavelength, the scattered 
field is predominantly in the forward direction (i.e. 

|/(fco|fco)| ^ l/(— ^o|^o)|)- Our forward approxi- 
mation consist in keeping only the contribution of 
the amplitude of diffusion /(fcjfco) in the direction 
of the incident wave kg . We write in using the hy- 



pothesis (|S5|) : 

h,oiKek\Keko) = 47r7(fc|feo) , (71) 

~47r7(fco|feo), (72) 

= 47r(T-fcofco)/(i^e,c^), (73) 

where 

f(K,,u;) = ^Sm = ^S2iO), (74) 



with 51(0) = 52(0) given by the Mie theory [1 H 
1^ "ss] . It is worth mentioning that the approxima- 
tion (|72|) is also valid for small scatterers (Rayleigh 
scatterers). In this case, the scattering amplitude 
/(fc|fco) doesn't depend on the direction of the in- 
cident and scattered wave vector k and fcoi since 
we have 

te,o{k\ko) ^ teA^)I ■ (75) 
From equation H68() . we show that: 

7(fc|fco) =7(feo|fco). (76) 

and we also obtain the coefhcient f{Ke,oj): 

ATTf{K„uj)^U^o{uj). (77) 

Furthermore, we see from equation (|52|l . that to 
zero order in density: 

Ce,o{k\ka)=te,o{k\ka), (78) 

and that the forward approximation H72|) applied 
to the operator Ce.o(fc|feo) in the low density limit. 
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We will suppose that the forward approximation is 
valid whatever the order in density for the operator 
Ce,o(fc|fco) E^nd we write: 

Ce,o(fc|fco) 2i Ce,o(feo|feo), (79) 

^(l^koko)CU\\ko\\,^). (80) 

With this hypothesis, only the path of type 1 in the 
figure are considered. This approximation also 
implied that the operator Ce,o (fc|fco) is transverse 
to the propagation direction feg. 




FIG. 1: Two difFerent paths which contribute to the mass 
operator S(||feo||,aj) = nCe,o{ko\ko) in the (QC-CPA) ap- 
proach. Only the path of kind 1 are taken into account in our 
forward scattering approximation 18U1 . 

From the previous hypothesis and the QC-CPA equations 
(|52|l . we obtain an equation on Ce,o (fco|fco): 



Ce^o(feolfco) = 47r/(/^e,c^)/±(fco) 
+4:Trn f{Ke,oj)rn{kQ) ■ Ce,o(fco|fco) 

where we have introduce the notation: 
I±{iio) = (-f - feofeo) , 

d^fcl 



m{ko) 



. ,„ h{ko~ki) ■I±iko)-G, (fci)-/i(fco) 



(81) 
(82) 
(83) 



Then, we have: 



Ce,o(fco|fco) = /±(feo) - 47rn/(A'e,t^)m(fco) 

■I^{ko)47rfiK,,Lo). (84) 

In using the classical properties of the Fourier transform, 
we write: 

m(fco) =7i(fco) .y dVe-"=« '-/i(r)Gr(r) •/±(fco) ■ 

(85) 

where we have used the translation invariance of the 
green function: Gg (r — rg) = (r,ro). We know 



that the Dyadic Green function has a singularity which 
can be separated in introducing the principal value of the 
Green function [ll[llil|13: 



GTir) = P.V.G,{r)- 



1 



■5{r)I. 



(86) 



where the principal value is defined by: 

P.V. I d3roGr(r-ro)-0(ro) 



= lima^o I d^roG^{r - ro) ■ 4>{ro) . 

'Sa{r) 



(87) 



with (p(ro) a test function and Sa{r) a spherical volume 
of radius a centered at r. This pr i ncip al value can be 
easily calculated, and we obtain hahflmol: 



P.V.Geir) 



47r||r|| 



1 



1 



iKe r 



3 



iXJIrll m\r\ 



In using polar coordinate in the integral (j85|l : 



d> . = 



47r 



(89) 



and the integral on solid angles given in Appendix ^ we 
obtain the following result: 



m(fco) 



m 

3i^2 



with 

m(ife) 
p{x) 



-t-oo 



rdrp{Ker) [g(r) — l]e 



(90) 



(91) 



smx 




sm X cos X 


X 


-{ 


x^ x^ 




1 




V IX 


x"^ 





sm X cos X 



(92) 

where we have assumed that ||A;o|| = K^. since the mass 
operator (|51|) is evaluated with this value in the equation 
(|60|l to obtain the effective permittivity. As the parti- 
cle cannot interpenetrate, we have (/(O) = and then 
/i(0) = -1. With equations ^fiOI25l51l84l90|l . we derive 
an expression for the effective wave number K^'- 



Kl = K\ 



(93) 



where the scalar m(Ke) is defined by equations (|91I92() 
and the scalar f(Ke,uj) is the forward scattering ampli- 
tude: /(fco|fco) — {I — koko)f{Ke,ijj) for a particle of 
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permittivity = Cg — ei + Ce within a medium of permit- 
tivity Ee • The relationship between the effective wave 
number Kf^ and the effective permittivity te is given by: 



ee{uj) = KllKl,,. (94) 
The formula H91I94|I are the main results of this paper. 

V. RAYLEIGH SCATTERERS 

We now show how to recover the low-frequency limit 
of the QC-CPA approach. First, we have to find an ex- 
pression for the T matrix for a single scatterer when its 
size is small compared to a wavelength (Kyacrs <Si 1). 
The T matrix te,o verifies the following equation: 

*e,o(T', To) = Ve,o{r) 6{r ~ro)+ J d^ri Ve,oM 

■G'^{r,ri)-te,o{ri,ro)- (95) 
where the potential is defined by: 



■ye,o(r) 



ei)Qs{r)I. 



(96) 



If we extract the singularity of the Green dyadic function 
Gg {r,ri) in using equation we obtain: 



*e,o(T", ro) = Vdipir) 6{r - Vq) + / d-^ri Vdip{r) 



where 



Vdzpir) 



[P.V.GT{r,ri)]-t,,o{ri,ro): 

Ve.o{r)' 



(97) 



3A'2 



with 



— ^vide Ctdip -* i 



47r 



^dip — 3 6e ^ 



£1 



(99) 

(100) 
(101) 



It's easy to recognize that the coefficient adip is the po- 
larization factor of a dipole. Hence, the singularity in the 
Green dyadic function describes the depolarization factor 
due to the induced field in the particles. The relationship 
between the singularity of the Green function and the de- 
polarization field acting on a particle has been described 
in numerous works [II [S Hi HI [H H HI HI • From 
the meaning of the coefficient adip, we inferred that equa- 
tion H97() describes the multiple scattering process by the 
dipoles inside the particle (where Qdir) ^ 0). 

As the particles are small compared to a wavelength, 
we use a point scatterer approximation: 



Qsir) 



6{r) 



(102) 



and the potential H99() becomes: 

Vd^p{r) ~ K^.de adipS{r)I. 



(103) 



In introducing the approximation H103() in equation l|97|) . 
the Dirac distribution provides an analytical answer for 
the T matrix of a single particle: 



*e,o(r", ro) =S{r - ro) 6{ro) te,o{i^) 



(104) 



te,o{^) = Kl^de adip I - Kl^de (^dip P-V-G^ (r = 0) 

(105) 

The principal value of the Green function at the origin 
can be evaluated in using a regularization procedure (60l 
El: 



P.V.G, (r = 0) = 



At 

6 TT 



6 TT 



(106) 



where the term At is proportional to the inverse of the 
real size of the scatterer [3, • Finally, the T matrix for 
a single particle is: 

t,^o{ri,r2)^5{r^-r^)5{r^)t,^o{K,,Lo)7, (107) 

^vide '^dip 



te,o{Ke,Uj) 



(108) 



It has been shown that the T matrix H107I108|) verifies the 
optical theorem and can present a resonant behavior due 
to the At term [3,|63|- The vahdity of the optical theorem 
is an important point, to insure that that the attenuation 
of the coherent wave due to scattering is correctly taken 
into account in the (QC-CPA) approach. Hence, the ex- 
pressions (|108|) must be used rather than the usual T 
matrix for a Rayleigh scatterer te^oiKe^co) = K'^ide'^dip 
which doesn't verify the optical theorem Further- 
more, from equations 168|) . we notice that: 



(109) 



The small size of the scatterers allow us also to approx- 
imate the term m{Ke) in equation H93|) . In fact, as 
there is no long range correlation in a random medium 
{g(r) — 1 ~ for r ^ Ts) and as if g ^ 1, we can 
evaluate the function p{Ker) in the integral H91|) for K^r 
close to zero. In using the limit: 



p{x) = - + o{x) X . 



(110) 



we obtain the following leading term of the real and imag- 
inary part of m{Ke)'. 



2 r+oo 

■miKe) = 3 y [g{r) - 1] 

r^dr[g{r)-l] + ... . (Ill) 



2iK, 

^ JO 
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If we keep only these two terms, equation ()93|l becomes 
in using the resuhs H1U9I1U8|) : 



(e, - ei) (1 - - Ke^' K^ac Ts)^ [f^ + i^«2]) + 3 ' 

(112) 

with fy — nvs the fractional volume occupied by the 
particles and wi, W2 defined by: 

r+oc 

47rn / rdr [g(r) - 1] = wi - At , (113) 
Jo 

Airn I r^dr[g{r)-l]=W2~l, (114) 



For non resonant Rayleigh scatterers, we can ne- 
glect the term At — 0, and also neglect the term 
^{el^"^ Kyacrs)'^ ^ compare to 1 — /„. Moreover, we 
have usually i?e(ee) ^ Im{ee) and equation pi2|l can 
be simplified into: 



ei + 



3(es - ei) Ee /, 



vol 



{(.s - ei) (1 - fvoi) + 3ee 



.2{KMerdY {^s - ^lY ^J"^ fvol 

[[es - ei)(l - f^oi) + 3eeJ^ 

where a Percus-Yevick correlation function for g(r) 
gives mill: 



W2 



(l + 2/„)2 



(116) 



The equation (|115|l is the usual low-frequency limit of 
the QC-CPA approach obtained by Tsang et al. [isl . 
In particular, we see that in the static-case (w = 0) the 
imaginary term in the right hand side of equation IjllSfl 
is null, and if we replace the effective permittivity Ce by 
ei in the right-hand side of equation H115|) . we recover 
the classical Maxwell Garnett formula: 



3(es — ei) ei f^oi 



{(.s - ei) (1 - fvoi) -I- 3ei 
which is usually written in the following form: 



Ee - ei _ - ei 

— Jv 



2ei 



2ei 



(117) 



(118) 



In comparing equation H117|l and (|112|l . we see that the 
scattering process modified the Maxwell Garnett formula 
by adding a new term: 



(119) 



whose imaginary part describes the attenuation of the 
coherent wave, and then, the transfer to the incoherent 
part due to the scattering of the wave. 



VI. KELLER FORMULA 

We are now going to show that the relation p3l) that we 
have obtained contain also the Keller formula |5^. This 
formula has recently been shown to be in good agreement 
with experimental results for particles larger than a wave- 
lengh 53, "5^. The Keller formula can be obtained in con- 
sidering the QC-CPA approach in the scalar case p^lT^ . 
The equations are formerly identical to equation H51l58fl 
where the dyadic Green function Gg : 



- VV 



\Kc ||r-ro|| 

47r||r|| 



(120) 



have to be replaced by the scalar Green function G° 
given by: 



Gr(r,ro,c^)- 



i -R'e 1 1 f—fa I 

47r||r|| 



(121) 



The first iteration of the scalar version of equation H52|) 
gives in using equation HA1() : 



Ki 



+ (471)2^2 fiK^Lu) 



+ 00 



sin Ke r 



(122) 



As was shown by Waterman et al, this development is 
valid if the following condition is verified: 



(4^)2n|/(X„c^)|2/i^e«l 



(123) 



In the geometric limit, the scattering cross section as for 
a single particle is in good approximation given by as — 
2 7Tr^. As the cross section is connected to scattering 
amplitude by the relation ag = ^\f{Ke,uj)\'^ and as the 
maximum density is n — 1/ws, we see that for particles 
larger than a wavelength the condition 1123|) is satisfied: 



{ATr)^n\fiK,,u)\yK, ~ 1/K,rs « 1 



(124) 



The equation H122|) . which has been derived by 
Keller 52], has proven to be in good agreement with 
experiments for particles larger than a wavelength. If 
we now use a Taylor development in equation H93|) . we 
obtain 



3Ki 



(125) 



This development is valid if the condition (|123|l is sat- 
isfied. In the geometric limits, we can approximate the 
function p(ii'e r) in the definition (|91|l of m{Ke), since for 
Ke r ^ 1 we have from equation 1)92(1 : 



p{x) 



smx 



a; > 1, 



(126) 
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The relation H125(l becomes: 



dr — ^ [5(7- 



l]e'^' 
(127) 



We see that equation (|127|l differ from the equation (|122|l . 
only by the factor 1/3K^, which is due to singularity 
of the vectorial Green function and consequently cannot 
be derived from the scalar theory developed by Keller. 
We also remarks that that to solve numerically our new 
equation l|9tj|l . we can use the same procedure that is used 
to solve the original Keller formula (|122() with the MuUer 
theory [s^. Consequently, we have derived a numerical 
tractable approximation to the (QC-CPA) scheme. 



the (QC-CPA) approach, we have added a far-field and a 
forward scattering approximations to (QC-CPA) scheme. 
In the low frequency limit, equation is identical with the 
usual result obtained under the (QC-CPA) scheme, and 
'in the high frequency limit the expression include the gen- 
eralization, in the vectorial case, of the result obtained 
by Keller. Further study is necessary to assess the lim- 
itation of this approach on the intermediate frequency 
regime. 



APPENDIX A: APPENDIXES 



sin 1 1 fen 



4-R 



1^0 



(Al) 



VII. CONCLUSION 

The intent of this paper has been to establish a new 
formula for the effective dielectric constant which charac- 
terize the coherent part of an electromagnetic wave prop- 
agating in a random medium. The starting point of our 
theory has been the quasicrystalline coherent potential 
approximation which takes into account the correlation 
between the particles. As the numerical calculation of 
the effective permittivity is still a difficult task under 



471 

cos||A;o|| Ih 



sin||fco||||r|| 
l|fcolPI|r||3 



+2 



llfcolPlklP _ 
cos||fco|| ||r| 

llfcolPIIHP 



(/ - fcofco) + 47r 
.sin||fco||||r|r 



sin||fco|| \\r\ 



\\ko\ 



(A2) 
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